Quasiclassical Surface of Section Perturbation Theory 
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Perturbation theory, the quasiclassical approximation and the quantum surface of section method 
are combined for the first time. This solves the long standing problem of quantizing the resonances 
and chaotic regions generically appearing in classical perturbation theory. The result is achieved by 
expanding the 'phase' of the wavefunction in powers of the square root of the small parameter. It 
gives explicit WKB-like wavefunctions even for systems which classically show hard chaos. We also 
find analytic solutions to some questions raised recently. 
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PACS: 05.45.+b, 03.65.Sq, 72.15.Rn 

Perturbation expansions in a small parameter e are im- 
portant in both classical and quantum physics. Not only 
are valuable approximations produced, but the break- 
down of the expansion can signal new physics. 

Poincarc found that classical perturbation theory [PT] 
on an integrable system fails in two [or more] dimensions 
for any e, due to 'smalLdenominators'. The Kolmogorov- 
Arnol'd-Moser [KAM]El theory greatly illuminated the 
subject and showed that the breakdown of PT signals 
chaos. Phase space trajectories of an integrable system 
lie on invariant tori. Under perturbation, periodic or- 
bits, on 'rational' tori, are destroyed except for one or 
more stable and unstable orbits. The rational tori are 
labelled by p, q, where p is the winding number and q is 
the number of returns to the SS per period. New invari- 
ant tori are formed around the stable orbits while chaos 
develops near the unstable orbit. The original tori near 
the rational one are also destroyed, to a width in action 
^feSpq. The characteristic action S pq generically drops off 
rapidly with q. This is usually pictured, as in Fig.l, on 
a surface of section [SS], a slice through the tori, where 
the structure of alternating stable and unstable orbits is 
called an 'island chain' or 'resonances'. 

Quantization of such a system has been of great inter- 
est. A rule of thumb is that only phase space structures 
of area Planck's constant h or greater are reflected in the 
quantum result. Thus, if \feS pq « Ti ordinary quan- 
tum perturbation theory works well. Above the 'Shuryak 
border' [SB]0, yfeS pq > h, ordinary perturbation theory 
breaks down as a number yfeS pq /Ti unperturbed quan- 
tum states are strongly mixed by the perturbation. Thus 
quantum perturbation theory for small h depends criti- 
cally on the relation between H, y/e, and the torus p, q. 

Small h is not a perturbation: rather the quasiclassi- 
cal approximation [QCA] is used. Combined QCA and 
PT has been studied aver the years: the perturbedrjrace 
formula most recently^. The SB does not appear inta, be- 
cause effectively only short times are considered, where 
classical perturbation theory does work. 

We here combine for the first time, PT, the. OCA, 
and the powerful QCA-SS method of BogomolnjoEj. We 
achieve quite complete and explicit results. Namely, we 
are able to find all the energy levels and wavefunctions in 
a WKB approximation for small e, provided h is not too 
small in a sense to be specified. Interestingly enough, 
we are led to express the wavefunctions in terms of a 



phase G which is generically expanded as a series in y/e 
rather than e. This series is generated in a novel recursive 
way: A partial n — l'th order solution is obtained which 
allows a partial n'th order solution which in turn gives 
the complete n — l'th order and a partial n + l'th order, 
and so on. We also give some numerical checks. 
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A number of problems suscep 
nique have appeared recentryu 
learned of several ongoing efforts 
related to a weakly deformed circular billiard. The 
Helmholtz equation (V 2 + A; 2 ) $ = is to be solved for 
eigenf unctions $ = f n and eigenvalues k = k a with, 
say, Dirichlet conditions ^f a (r, 0) = on the bound- 
ary, dB. The latter is expressed in polar coordinates by 
riff) = Rq + eA.R(9). These and similar boundary pertur- 
bations have heretofore been treatedliJ by methods valid 
only below the SB. 

We illustrate with the stadium billiardB which has two 
semicircular endcaps of radius Rb connected by paral- 
lel straight sides of length 2a. Then e = a/Rs is as- 
sumed small and Rq is taken as Rq w Rb + 2a / V while 
AR/R w | sin 1 - 2/tt. This 'stadium' choice of AR has 
a discontinuous first derivative so KAM does not apply. 
Fig. la shows that the classical map deviates from the 
new invariant tori [given approximately by Iwkb(&) de- 
fined below] after about 1/y/e iterations. We also show 
results, Fig. lb, for a 'smoothed stadium', a truncated 
Fourier series of the 'stadium' AR(9), where the orbit 
stays on a new invariant torus. 

Classically the stadium is chaotic with no stable ac.- 
bits. Orbits diffuse in angular momentum at long timesu. 
It was thought that such hard chaos systems do not 
have simple, analytically expressible wavefunctions when 
quantized. Thus, qualitative and statistical questions, 
such as the existence and statistics of localization, have 
been considered. Our explicit analytic results were there- 
fore quite unexpected, and we are able to interpret the 
results directly in terms of analytic wave functions. The 
results are possible because it is the short time, nearly 
regular behavior which determines the quantization. 

In quantum language, we take units Rq = 1, ft = 1, 
particle mass = 1/2, so k is the dimensionless wavenum- 
ber, [equivalent to 1/K\. We take the billiard boundary 
dB as SS. Then Bogomolny's unitary operator iaa 
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T(9,9',k) 



k d 2 L(0,0') 
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exp (ikL (0, 0')) (1) 



where L is the chord distance between two points, spec- 
ified by polar angles, on 8B. Expanding, 



kL(9,9') = 2k 



1 + e 



AR{0) + AR{0') 



= k(L + eL 2 + ....) (2) 

The energy levels [k a ] of the system are given in QCA§ 
by solution of dct(l — T(k)) = 0. Our seemingly more 
difficult technique studies 



i>{6) = / d9'T(9,9',k)i/}(9') 



(3) 



solvable only for k — k a . [ijj ps d^/dn on 8B] 

We start with the Ansatz ip(0) = exp (iaf(9)) where 
df/d6 = f ~ 1 and k » a » ke. This Ansatz 
represents a superposition of angular momentum states 
\l\ ~ a << Z max = fc, and for a > 1 conveniently ex- 
presses the mix of states needed to diagonalize the Hamil- 
tonian above the SB. 
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FIG. 1. SS [angular momentum vs angle] of orbits for a) 
'stadium' and b) 'smoothed stadium' e = 0.01. Points on 
Iwkb{9), for 'continuum', O, 'separatrix',+, and 'bound', *, 
values of E m , respectively. Three orbits, each iterated 1000 
times, coalesce into solid lines in b), where KAM applies. 
Orbits started at the symbols iterated forward and backward 
15 times appear as dots in a) where KAM fails. Only short 
time structure is regular, c) Husimi plots for exact states of 
Fig. 2, as well as a 'scar' state '032. Square has area h. 



According to the stationary phase [S$] method, the 
0' integral is dominated by the region 0' ~ ± ir where 
kL — 2k | sin \{9 — 9')\ is stationary. Expand sin |(0 ~ 
9') rs 1 - h89 2 , 69 = 9' - 9 - tt to find 



S(9,9') = kL « 2k- 



\k59 2 



ke(AR{0) + AR{9')) (4) 



[In Eq.(||) we replaced Lq by 2, its stationary value, when 
multiplied by e.] 

Regarding S as a classical generating function, we ob- 
tain the surface of section maps (/', 0') — ► (/, 0) found 



earlier^ by I = QS/89, V = -8S/80'. Motivated 
by this, Borgonovfl has studied the T operator and 
classical map given by Eq.(Q) with AR = \sm9\ and 
59 = 0' — . This system is 'almost' the quantum 
kicked rotor-classical standard mapllil, which corresponds 
to AR = sin (9. Thus, in addition to solving the distorted 
billiard problem, we can also solve an important class of 
quantized perturbed twist maps. 

Returning to Eq.(p3), we expand all functions of 0' 
about 9 + ir. I.e. /(#') re f(9 + tt) + 59 f'(9 + n) to 
order 59, [since a « k] and AR(9') pa AR(9 + tt), since 
a >> ke. Doing the integral reduces Eq.(^) to 

exp[ia/(0)] 

= % exp [i (2k + [aff/k + keV(9) + af(9 + tt))] (5) 

where V(9) = AR(9) + AR{9 + tt). 

For Eq.rtq) to hold, the exponents of order a must com- 
bine to give a constant c, i.e. f(9 + n) = f (9) + c. 
Now take a — ky/e so a solution is possible provided 
(f) 2 + V(0) is a constant, which we call E rn . Thus 



f(9)=± / d9'^E m - V{9') 
Jo 



(6) 



reminiscent of elementary WKB theory. The constant of 
integration is irrelevant. Notice V(9) = V(9 + tt) => 
fiS + tt) = / (0) + c. We define I WK b(9) = kyfef(9). 

Assuming for now that E m > V , [a 'continuum' state], 
we must choose E m such that kb J Q * d0' \J E m — V(0') = 
2Txm where m is integer, b = i/e and soc = Txmjkb. The 
condition giving the energy is 

exp [i (2k + kb 2 E m + kbc + tt/2)] = 1 = exp(27rin) (7) 



. For AR = 0, this reduces 
-j)27r equivalent to Debye's 



which has solutions k = k 
to 2k + m 2 jk + nm = (n 
approximation to Bessel's function, valid for k large and 
m/k small. Thus, this Ansatz produces states labelled 
m, n with m an integer angular quantum number satis- 
fying |m| << k n>m ps 7T7i. There are three symmetries, 
reflections about the two principal axes and time rever- 
sal, which allows real wavefunctions. Thus the even-even 
states are ip = cosaf(9) and m is an even integer. [We 
choose the lower limit in Eq.(|^) to be at a minimum of 
V.] This result allows an explicit estimate of tpi [angular 
momentum representation] which, for |Z| > k^/e, decays 
exponentially for smooth V and as l~ A for the stadium 
case. This localization was firstQ thought to be dynami- 
cal localization analogous to Anderson localkationE3, but 
nowQ [for ke 2 < 1] is attributed to CantoritH 

If E m — V changes sign there are 'bound state' regions 
near the minima of V [at 9 = 0] where E m > V, e.g. let 
the region be \9\ < m < tt/2, where = m is a 'classical 
turning point' of the motion. The even-even quantization 
condition is now, approximately, cos(af(9 m )) = 0, or 
af(9 m ) = (m + !)tt and tp ps 0, 9 m < \9\ < tt/2. In 
this approximation there is a degeneracy between even- 
even and even-odd symmetry. This treatment neglects 
tunnelling into the forbidden region V > E m , as well as 
effects on the amplitude of the wavefunctions. 
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The bound states quantize the stable resonance islands 
and the continuum states quantize the unstable and per- 
turbed KAM regions. A minimum in V is at a stable pe- 
riodic orbit, and a maximum at an unstable one. More 
correctly, if V does not have sufficiently many deriva- 
tives, the stable orbits can disappear, but the quantum 
system is hardly affected, if the wavelength is not too 
short. 'Scars' of unstable orbits, Fig. lc, are states with 
E m just greater than the maximum V. Fig. 2 shows a 
WKB state and two indistinguishable numerically ob- 
tained states, all with the same value of fc^/e, one with 
ke = 1.8 the other for ke = 0.18. This shows the state 
depends dominantly on k\fe and suggests no qualitative 
changes occur at ke » 1. Husimi plots of these states are 
shown in Fig. lc. 

Borgonovii has numerically calculated an average lo- 
calization width in angular momentum, l a , where, [in 

effect] I 2 = J^a I Co. 1 2 J d9 \ip' a \ 2 and c a is the normalized 
zero'th Fourier component of the eigenstate ip a . Accord- 
ing to the results just obtained, the c a 's should be rel- 
atively small for 'continuum' states, [Fig. 2] since the 
phase is not stationary, so the 'bound' states dominate. 
Then ll » a 2 J2a\ca\ 2 Jdd(E m{a )~V(d)) \smaf a (0)\ 2 
The sum is now a sort of average (E — V) which is of 
order unity and nearly independent of a. Thus l a cx k^/e. 
Borgonovi fixed k and increased e, agreeing with this re- 
sult until ke 2 ps 1 . We show below that our theory should 
fail at that point. We note that the result depends on 
the definition of the c a 's. The result can be quite differ- 
ent if the c a 's are chosen to be the overlap of ip a with 
some high angular momentum state, for example. Fig. 
2 shows a high angular momentum state, away from a 
resonant torus, which has a much narrower distribution. 
[See the next paragraph.] 
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FIG. 2. 'Stadium' potential |sin# vs angle. States and 
potential are symmetric about zero angle. Bound and contin- 
uum, WKB and exact states are shown, with zeroes at WKB 
'energy' parameter E m . fee 1 / 2 = 42.3 is fixed. Inset: Angular 
momentum representation of continuum state m = 48 and 
exact state near angular momentum m = 168. 



We turn to general angular momenta and higher orders 
in e. We look for solutions of the form i/j — exp(iG(0')), 
where G = W + k(ef 2 + e 2 f A + ■■■)■ The /'s are 2tt- 
pcriodic and I < k is integer. This, if successful, is a 
usual PT for G. The S<$> angle is 0' w + 8j where 
0; = — 2sign(7) cos -1 (7/fc). Expanding as before the order 
ke condition is 

f2(e + Qi)-f 2 (6)=L 2 (9). (8) 

We use l 2 and L 2 {0) as the constant and variable parts of 
L 2 (9,9 + Qi). The constant part kel 2 contributes to the 
phase of Eq.(Q). Eq.(||) is solved in terms of r ^ Fourier 
components, i.e. f 2r = (exp(ir0;) — 1) 1 L 2r . This a 
good solution unless the denominator is excessively small. 
It never strictly vanishes since Qi/2tt is an irrational 
number. However, if 0; is close to Q pq = 2?rp/q, where 
p/q is a rational number, [corresponding to the strongly 
perturbed rational tori of classical perturbation theory] , 
then the denominator will be small if r is a multiple of q. 
It will still be a good solution if L 2r vanishes or is suffi- 
ciently small. Generically L 2r decreases rapidly for large 
r. If the small denominators are thus compensated by 
small numerators, this perturbation theory can be car- 
ried to higher orders by the methods described below. 
If not, we need to refine the approach along the lines of 
our first Ansatz which corresponds to q — 2. This small 
denominator problem is the analog ia, QCA of the small 
denominator problem of classical PTEJ. 
We are thus motivated to consider 

i> = exp [i (l pq d' + k{bh + b 2 f 2 + b 3 f 3 + ...))] (9) 

The [non-integer] angular momentum l pq is chosen to 
make the stationary point 6' = 9 + po . Expanding as 
before, the order b requirement is fi(9 + <d pq ) — /i(0) + 
IpqOpq/kb = c =constant implying /{ is g-periodic, i.e. 
periodic with period Q pq . At order b 2 we have 

S;, 1 {fif + L 2 {0) + h(9 + O pq ) - f 2 (6) = E m (10) 

where L 2 (9) is the variable part of L 2 (9, 9 + Q pq ), S pq = 
|sini0 po | and E m is to be determined. We divide 
Eq.(|l0[) into q-periodic and non q-periodic parts. The 
nonperiodic terms f 2 and L 2 must combine to give a q- 
periodic result, thus 

f2(9 + e pq )-f 2 (9)+L 2 (9) = V q (9) (11) 

where V q (9) is to be determined. We 'g-average' both 
sides giving V q {9) = | J2]=i L ^{9 + j@ Pq )- Expressed in 
Fourier components, V q (9) = '^Zi L 2q ie lql6 and f 2 (9) = 
Yli ( 1 - e im " q ) ~ 1 L 2l e lW + f 2 (9) . The prime indicates that 
integers I divisible by q are not included in the sum and 
f 2 (9) is an g-periodic function not yet determined. Then 

/x(0) = ±SH 2 J d9'^E m - %{9>) (12) 

and considerations like those discussed earlier for q = 2 
fix the quantization of E m . The size of V qi which de- 
creases rapidly with q, determines if powers of y/e~ rather 
than e are needed. 
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Order b 3 is more complicated: L2{0,9') and fa are 
expanded to 59, fa to S9 2 and Lq to S9 3 . The integral of 
Eq.(|) is thus 



d69 exp 



— ife 
~3T 



r" 



ik 



+ -(L^ + b.fa)59 2 +tF'89 



(13) 



where F' = fcb/f + fc6 2 F 2 ' with F' 2 = f 2 + L' 2 . We denote 
derivatives evaluated at 9' = 9 + <d pq by primes. [This 
integral is done over a region near the original stationary 
point. The new stationary point coming from S9 3 is not 
meaningful.] The width of contributing angles 69 is of 
order fc -1 / 2 which is small. However, the shift of the 
center of the contributing region is expressed by a power 
series in b whose leading term is — 6/{/Lq. If kb 2 > 1, the 
shift cannot be neglected. Thus, to order b 3 we require 
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where C3 is a constant. Let F 2 = fa 2 + A(9), where A has 
already been determined by lower order considerations. 
Eq.(pT[) can only be satisfied if the g-average of the left 
hand side vanishes. This determines f 2 by 
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This expression must also have vanishing angular aver- 
age, since f 2 is the derivative of a periodic function, which 
determines C3. Thus fa is determined up to an irrelevant 
integration constant, and, then as before, fa is deter- 
mined up to an g-periodic function. 

If kb 3 << 1, we may stop here. If not, we can 
continue finding higher order corrections, expanding to 
higher powers of 39 and keeping the terms L 4 , L 6 , ... in 
the expansion of the phase of the T operator. The series 
will be effectively terminated at order n when kb n << 1. 
However, the method may break down sooner, indicating 
a change in the fundamental physics. 

For example, 'bound state' solutions of Eq.(|l2|) give 
infinite second derivatives /{' when the square root van- 
ishes, i.e. at the 'classical turning points'. This, however, 
can be taken into account to give the familiar turning 
point corrections of elementary WKB theory. 

In the case AR = |sin 6*| , there are 6- function singular- 
ities in /{" and L 2 . These large derivatives invalidate the 
expansion. Thus in the Bunimovich problem we expect 
our solution to break down when ke 2 > 1. Numerical 
resultsu seem to confirm this expectation, giving two dif- 
ferent behaviors in either side of this border. 

In principle, we can use this technique to study pertur- 
bations of any two dimensional integrable system. 'Sim- 
ply' use action angle coordinates fa, I2, Oi, ©2, and take 
as surface of section ©1 = 0. The T operator will have a 
phase fc(S , o(6 2 -6 2 ) +£5*2(62,02) + •■■) and the rest is 
pretty much the same as above. Other coordinates may 



be more convenient in practice, however. The circle is 
nice because the action-angle coordinates are immediate. 

There are other applications of this technique in non- 
perturbative settings, in which certain classes of eigen- 
states can be found. The germ of the method-first ap- 
peared in the study of the ray splitting billiardtfl, and it 
can be used to find the well known 'bouncing ball' states 
in the [large e] stadium billiard. 

We have thus produced a fairly general theory allowing 
us to find the effect of perturbations on integrable quan- 
tum systems which exploits the quasiclassical approxi- 
mation and the surface of section technique. If the per- 
turbation classically gives rise to resonances big enough 
to influence the quantum problem, we must expand in 
the square root of the small parameter. If the resonances 
are small, a simpler expansion works. 
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